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Analyseren van opeenvolgende datasets 



De uitvinding heeft betrekking op een werkwijze vooi het analyseren van 
opeenvolgende datasets waarin de afzonderlijke datasets dataelementen omvatten die 
datawaarden toevoegen aan ruizntelijke posities, 

5 Een dergelijke werkwijze is bekend uit het Amerikaanse octrooischtift 

US 5 431 161. 

De bekende werkwijze maakt gebruik van datasets die zijn verkregen met een 
Single Photon emission comuted-tomography ASPECT) afbeeldingsmetbode. Deze bekende 

10 werkwijze is vooral bedoeld voor diagnosticeren van hartaandoeningen, met name van 
ziekten van de kransslagaderen. Met de SPECT afbeeldingsmetbode worden SPECT- 
afbeeldingen van doorsneden van het hart van een te onderzoeken patient gemaakt 
Vervolgens dient de gebruiker van de bekende methode uit de SPECT-afbeeldingea 
handmatig gebieden van belang te selecteren. De gebieden van belang zijn met name de 

1 5 delen in de SPECT-afbeeldingen die betrekking hebben op het myocardium (de hartspier). 
De bekende werkwijze bepaalt vervolgens de perfusie door het myocardium om een indruk 
van de mate van doorbloeding van de hartspier te krijgen. Hiertoe wordt de verhonding van 
de gemiddelde helderheden van de helderste gedeelte van SPECT-beelden in de eind- 
systolische en de eind-diastolische fase van het te onderzoeken hart bepaald. 

20 

Een doel van de uitvinding is een werkwijze te verschaffen waarmee een 
nauwkeuriger analyse van de datasets mogelijk is. In het bijzonder is een doel van de 
uitvinding een werkwijze te verschaffen waarmee het gebied van belang in de opeenvolgende 
datasets preciezer kan worden opgespoord. 

25 

Dit doel wordt bereikt met een werkwijze voor het analyseren van 
opeenvolgende datasets waarin 

- een locaal intensiteitsverloop [/fo#]woidi opgemaakt uit datawaarden in opeenvolgende 
datasets in overeenkomstige raimtelijke posities en 
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_ ^t een of meer van de opeenvolgende datasets een gebied van belang wordt opgespoord 

aan de hand van het locale intensiteitsverloop 
_ waarbij in net gebied van belang het locale intensiteitsverloop voldoet aan een vooraf 

gegeven eigenschap. 



Volgens de uitvinding wordt en in de opeenvolgende datasets het gebied van 
belang opgespoord aan de hand van het locale intensiteitsverloop van de datawaarden in de 
datasets. De uitvinding is met name goed bruikbaar wanneer elk van de datasets een 

10 afbeelding van een te onderzoeken object vertegenwoordigt. In dit geval zijn de datawaarden 
helderheidswaarden of grijswaarden van de pixels van de afbeelding. De uitvinding is 
overigens ook bruikbaar wanneer de datawaarden de plaatselijke waarde van een fysische 
grootheid vertegenwoordigen, hierbij kanhet gaanom temper atuur, druk, concentraties van 
bepaalde stoffen en dergelijke. Het intensiteitsverloop vertegenwoordigt de variaties van 

15 datawaarden in opeenvolgende datasets in overeenkomstige posities. Bijvoorbeeld, wanneer 
de datasets in de tijd elkaar opvolgen vertegenwoordigt het intensiteitsverloop de variatie in 
de tijd van de datawaarden in overeenkomstige posities. De uitvinding is bijvoorbeeld goed 
bruikbaar voor magnetische-resonantie perfusie afbeeldmgsmethoden; daarbij hebben de 
opeenvolgende datasets betrekking op opeenvolgende fasen van een perfusieproces. De 

20 datasets worden hierbij verkregen door reconstructie uit magnetische-resonantiesignalen die 
in de opeenvolgende fasen van het perfusieproces met een perfusiegevoelige magnetische- 
resonantie afbeeldingsmethode zijn opgenomen. De aanpak volgens de uitvinding berust op 
het inzicht dat er een verschil bestaat tussen het verloop in de tijd van de datawaarden in het 
gebied van belang en van de datawaarden buiten het gebied van belang. Zelfs komthet vaak 

25 voor dat er in het gebied van belang een bepaald karakteristiek intensiteitsverloop 

plaatsvindt. Zo'n karakteristiek intensiteitsverloop blijkt gemakkelijk op te sporen te zijn, 
met name door het actuele intensiteitsverloop te vergelijken met een vooraf gegeven 
intensiteitsverloop. Meestal is het voldoende deze vergelijking te beperken tot een of enkele 
karakteristieke eigenschappen van het intensiteitsverloop. 

30 Deze en andere aspecten van de uitvinding worden verder uitgewerkt aan de 

hand van de volgende uhvoeringsvoorbeelden die zijn gedefinieerd in de afhankeUjke 
Conclusies. 

Bij voorkeur wordt het locale intensiteitsverloop opgemaakt voor 
afzonderlijke blokken die betrekking hebben op grotere of ldeinere ruimtelijk gebieden, maar 
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bij voorfeeur meerdere datawaarden van de individuele datasets omvatten. Liefet wordt deze 
maatregel toegepast in situaties waarin het locale intensiteitsverloop aanmerkelijk verandert 
over aanmerkelijke niimtelijke afstanden, maar slechts gering varieert over kortere afstanden. 
Hierdoor is het niet nodig om het locale intensiteitsverloop voor elk van de afzonderlijke 
5 ruirntelijke posities te bepalen, zodat de hoeveelheid berekeningen om tot het locale 

intensiteitsverloop te komen wordt beperkt. De grootte van de blokken is liefst ongeveer de 
afstand waarover en aanmerkelijke verschillen in het locale intensiteitsverloop optreden. 
Goede resultaten wordt bereikt met niet-overlappende blokken van ongeveer 32x32 
beeldelementen (pixels). In deze afzonderlijke blokken wordt het plaatselijk gemiddeld locale 

1 0 intensiteitsverloop berekend. Hiertoe worden er in de afeonderlijke blokken telkens in de 
individuele datasets locale gemiddelden van de datawaarden berekend. Uit deze locaal 
gemiddelde datawaarden wordt het locaal gemiddelde intensiteitsverloop afgeleid De locale 
gemiddelde datawaarden worden per blok in de individuele datasets verkregen door filteren 
met een bij voorbeeld Gaussisch profiel. Hierdoor wordt vermeden dat toevallige fluctuates, 

1 5 zoals door mis, het opsporen van het gebied van belang nadelig beKhvloeden. 

Een geschikte karakteiistieke eigenschap van het locale intensiteitsverloop is 
de tijdafgeleide van het locale intensiteitverloop. In veel situaties is in het gebied van belang 
de tijdafgeleide van het locale intensiteitsverloop duidelijk verschillend van de tijdsafgeleide 
van het locale intensiteitsverloop buiten het gebied van belang. Met name is de tijdafgeleide 

20 van het locale intensiteitsverloop erg geschikt om toe te pas sen in de cardiologie. Met allerlei 
radiologische technieken, zoals magnetische resonantie angtografie, r5ntgenangiografie of 
computer tomografische angiografie is het mogelijk om opeenvolgende datasets op te nemen 
waarin het hart van een te onderzoeken patient wordt afgebeeld. De datawaarden in de 
datasets representeren dan locale intensiteitswaarden in de angiografische afbeelding. Zo'n 

25 afbeelding betreft bij voorbeeld een dwarsdoorsnede door het hart van de te onderzoeken 

patient, maar het is ook mogelijk een ruirntelijke (drie-dimensionale) afbeelding van het hart 
van de te onderzoeken patiSnt te maken. In deze toepassing zijn gebieden van belang 
bij voorbeeld de gedeelten in de datasets die betrekking hebben op het linker en zechter 
ventrikei en het myocardium van het hart van de te onderzoeken patient. Bij het opnemen van 

30 de datasets is er verder een contrastmiddel ^an de patient toegediend. Het locale 

intensiteitsverloop vertegenwoordigt dan het passeren van bloed met contrastmiddel door de 
afzonderlijke anatomische delen van het hart van de te onderzoeken patiSnt. Verder zijn er 
aanzienlijke verschillen tussen het locale intensiteitsverloop in de afzonderlijke anatomische 
delen van het hart. Bijvoorbeeld treedt eerst een aanzlenlijke snelle toename van de intensiteit 
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op gevolgd door een wat langzamere afiiame van de intensiteit op in het linker ventrikel. In 
het rechter ventrikel treedt een dergelijk intensiteitsprofiel op met een zekere tijdvertraging 
en nog later doet dat intensiteitsprofiel, maar met minder sterke tijdvariaties van de intensiteit 
zich in het myocardium voor. Zo lukt het volgens de uitvinding om nauwkeurig en 
5 automatisch uit een angiografische (cardio) dataset de linker- en rechterventrikels en het 
myocardium op te sporen. 

Met name doen zich aanzienlijke verschillen voor tussen de tijdafgeleiden van 
het locale intensiteitsprofiel in respectievelijke het linker- en rechter ventrikel, het 
myocardium en het gebied buiten het hart van de te onderzoeken patient Hierdoor is de 
1 0 tijdsafgeleide van het locale intensiteitsverloop erg geschikt om nauwkeurig en automatisch 
uit een angiografische (cardio) dataset de linker- en rechterventrikels en het myocardium op 
te sporen. 

In het bijzonder lukt het om het gebied van belang op te sporen aan de hand van blokken van 
d&taelementen waar de vaiiaties van het locale intensiteitsverloop, bij voorbeeld de 

1 5 tijdsafgeleide van het locale intensiteitsverloop een vooraf bepaalde plafondwaaide 
overschrijden en/of kleiner zijn dan een vooraf bepaalde bodemwaarde. Juist in de 
cardiologische tpepassing van de angiografie lukt het om zo bij voorbeeld het linker- en 
rechter ventrikel op te sporen. In het linker- en rechterventrikel treedt een snelle aanzienlijke 
toename gevolgd door een aanzienlijke afiiame van de locale intensiteit op, dat betekent dat 

20 er een overschrijding van de plafondwaarde door de tijdsafgeleide van het locale 

intensiteitsverloop gevolgd door een waarde van de tijdsafgeleide lager dan de bodemwaarde. 
Deze bodemwaarde is bijvoorkeur negatief. Zo lukt het om gedeelten in de datasets op te 
sporen waarin er een afhemende intensiteitsverloop optreedt. Zo'n afhemend 
intensiteitsverloop doet zich bij voorbeeld voor doordat de concentrate contrastmiddel in het 

25 bloed in het hart van de te onderzoeken patient daalt na de passage van een contrastbolus. In 
de pralctijk is het vaak voldoende om alleen het linker- en rechter ventrikel uit de 
opeenvolgende datasets op te sporen aan de hand van de tijdsafgeleide van het locale 
intensiteitsverloop. Het myocardium kan vervolgens aan de hand van de inmiddels bepaalde 
ligging van het linker- en rechterventrikel worden opgespoorcL 

30 Volgens de uitvinding lukt het om aan de hand van het locale 

intensiteitsverloop de dataelementen te classificeren als al dan niet behorende tot het gebied 
van belang. Bijvoorbeeld lukt het goed om aan de hand van het optreden van waarden van de 
tijdsafgeleide boven de plafond waarden en negatieve waarden van de tijdsafgeleide onder de 
bodemwaarde dataelementen te classificeren als behorend tot het linker- of rechter ventrikel. 
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In een verdere voordelige implementatie van de werkwijze volgen? de 
uitvin ding wordt het locale intensiteitsverloop vergeleken met een referentie 
intensiteitsverloop voor het gebied van belang. Zo'n gemiddeld intensiteitsverloop wordt 
bijvoorbeeld vooraf empirisch bepaald. Bijvoorbeeld worden empirisch van een groot aantal 
5 patiSnten de intensiteitsverlopen in de linker- en rechterventrikels en het myocardium 

bepaald. Van deze empirisch bepaalde intensiteitsverlopen voor de afzonderlijke delen van 
het hart worden individuele gemiddelde intensiteitsverlopen gevonnd. Deze gemiddelde 
intensiteitsverlopen zijn goed bruikbaar als referentie. 

De vergelijking van het locale intensiteitsverloop met het referentie 

10 intensiteitsverloop blijkt nauwkeurig uitgevoerd te kunnen worden door de tijdsafgeleide van 
het locale intensiteitsverloop te vergelijken met de tijdsafgeleide van het referentie 
intensiteitsverloop. 

De classificatie van dataelementen is nauwkeurig te doen zonder veel 
rekenwerk door het locale intensiteitsverloop voor het betreffende dataelement met het 

1 5 referentie intensiteitsverloop te vergelijken. Bijvoorbeeld wordt de vergelijking gedaan door 
de correlate tussen het locale intensiteitsverloop en het referentie intensiteitsverloop te 
bexekenen. Deze vergelijking kan vooral nauwkeurig en kwantitatief worden uitgevoerd met 
behulp van een gelijkheidsmaat Zo'n gelijkheidsmaat duidt de mate van overeenstemming 
van het locale intensiteitsverloop met het referentie intensiteitsverloop aan- 

20 In een voordelige implementatie van de werkwijze volgens de uitvinding 

worden sterke ruimtelijke gradienten van datawaarden in gedeehen van de datasets 
gemaskeerd. Met maskeren wordt in dit verband bedoeid dat zulke sterke ruimtelijke 
gradifinten worden weggelaten of voor de verdere bewerkingen genegeerd. Door zulke sterke 
ruimtelijke gradienten samen met gevolgen van beweging, bijvoorbeeld tijdens het opnemen 

25 van de datawaarden, kunnen er aanzienlijke kunstmatige temporele gradienten optreden. 
Door gedeelten in de afzonderlijke datasets met sterke ruimtelijke gradienten van de 
datawaarden te maskeren wordt nadelige invloed van zulke kunstmatige temporele gradienten 
op het opsporen van het gebied van belang vermeden. Het maskeren wordt bijvoorbeeld 
uitgevoerd door in de datasets plaatselijk gemiddelde gradienten te berekenen en de 

30 plaatselijk gemiddelde gradienten af tetrekken van de datawaarden in de datasets- Qnder de 
omstandigheid dat tevoren bekend is dat het gebied van belang zelf (ook) wordt gekenmerkt 
door sterke ruimtelijke gradiSnten kan daar rekening mee gehouden worden. Bijvoorbeeld bij 
MR-angiografie met een perfusie wordt het maskeren beperkt tot de datasets vooraf aan het 
optreden van perfusie en/of na afloop van de perfusie. 
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Liefst worden de datasets onderling gecorrigeerd voor verschijnselen die door 
beweging, bijvoorbeeld gedurende de opname van de datasets, ontstaan. Zo'n coirectie wordt 
wel als 'onderlinge registratie' aangeduid. In opeenvolgende datasets worden ruimtelijke 
veiplaatsingen uitgevoerd, dat wil zeggen de ruimtelijke posities waaraan de datawaarden 
5 zijn toegevoegd worden aangepast 20 dat overeenkomstige gedeelten in opeenvolgende 
datasets telkens bij dezelfde ruimtelijke posities liggen. Zulke overeenkomstige gedeelten 
hebben bijvoorbeeld telkens betrekking op eenzelfde (gedeelte van) een object dat door de 
datasets op achtereenvolgende momenten in de tijd wordt afgebeeld. 

De werkwijze volgens de uitvinding is verder geschikt om meerdere 

1 0 verschillende gebieden van belang nit de opeenvolgende datasets op te sporen. In het 
bijzonder lukt het goed om verschillende gebieden van belang op te sporen waarbij er 
verschiUen bestaan tussen de locale intensiteitsverlopen in deze verschillende gebieden. Dit 
doet zich met name voor in de cardiologische angiografie waarbij de opeenvolgende datasets 
betrekking hebben op de passage van bloed met contrastmiddel door de linker- en 

1 5 rechterventrikels en het myocardium van het hart van de te onderzoeken patient Deze 
afzonderlijke gebieden van belang worden nauwkeurig met de werkwijze volgens de 
uitv inding opgespooid. Bij voorkeur worden de overeenkomstige gebieden van belang in de 
individuele datasets onderling geregistreerd om kunstmatige verschijnselen te corrigeren die 
ontstaan door bewegingen, zoals door het kloppen van het hart en de ademhaling van de te 

20 onderzoeken patignt Volgens een voordelige implementatie van de werkwijze volgens de 
uitvinding worden er maximum intensiteitsprojecties in een volgoiderichting door de 
opeenvolgende datasets gemaakt van de afeonderlijke ruimtelijke gebieden van belang die al 
zijn opgespooid. Verschillen tussen deze afzonderlijke maximum intensiteitsprojecties 
vonnen thema-afbeeldingen waarin telkens een van de gebieden van belang versterkt wordt 

25 weergegeven. Met name lukt het om van afzonderlijke gebieden van belang versterkt ten 

opzichte van andere opgespooide gebieden van belang weer te geven in situaties waarin in de 
afzonderlijke gebieden van belang telkens een vergelijkbaar intensiteitsverloop, maar met een 
onderlinge veitraging in de tijd optreedt. Zo'n situatte doet zich met name in de 
cardiologische angiografie door waar perfusie van bloed met contrast middel eerst door het 

30 rechter ventrikel, dan door bet linker ventrikel en uiteindelijk door het myocardium passeert 
Door een ma ximum intensteitsprojectie door de opeenvolgende datasets van het linker en 
rechterventrikel te maken in de tijdspanne tussen de passage door het linker en rechter 
ventrikel te maken van de gebieden die op het linker- en rechter ventrikel te maken en af te 
trekken van een andere maximum intensiteitsprojectie op een moment na de passage door het 
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linker ventrikel door de opeenvolgende datasets van het gebied van belang dat betrekking 
heeft op het linker ventrikel oatstaat een thema-afbeelding waaiin het linker ventrikel en het 
myocardium versterkt zijn weergeven. 

Het blijkt dat de globale vorm van het gebied van belang vaak, vooral in cardiologische 
5 toepassingen wel vooraf bekend is, bijvoorbeeld is het gebied van belang min ofmeer rond. 
Uit de thema-afbeeldin g is met gebruik van de voorafbekende globale vorm de begrenzing 
gemakkelijk af te leiden. Met name is hiervoor interpolatie met een restrictie op de 
voorafbekende vorm geschikL Ook is interpolatie aan de hand van een zogenoemd statistical 
shape model geschikt am de begrenzing vast te stellen. Zo*n statistical shape model is op 

1 0 zich bekend uit TJ. Cooten et al. "Active shape models - their training and application", 
Comp. Vis- Image Underst61(l):38:59, 1995. en ook uit AJD. Brett and CJ. Taylor, 
"Automated construction of 3D shape models using harmonic maps f \ in Medical Image 
Understanding and Analysis, London, July 2000, pages 175-178. 

Door de thema-afbeelding waaiin het ronde gebied van belang versterkt is 

1 5 weergegeven naar poolcoSrdinaten te transformeren is de rand van het gebied van belang 
gemakkelijk op te sporen. Bij voorkeur wordt de oorsprong van het poolcoordinaten stelsel 
gekozen in het midden van het ronde gebied in de thema-afbeelding- In de naar 
poolcoordinaten getransformeerde thema-afbeelding treedt de rand van het gebied van belang 
op als een duidelijke overgang in de dalawaarden, die zich in een weergave uit als een rechte 

20 rand. Zo'n rechte rand is gemakkelijk op te sporen. In het bijzonder is het op deze manier 
gemakkelijk om in de thema-afbeelding het gebied van het myocardium nauwkeurig op te 
sporen. Vervolgens kan de plaatselijke doorbioeding van het myocardium bestudeerd 
worden. 

In een voordelige implementatie van de uitvinding wordt voor den of meer van 
25 de individuele datasets aan de hanrj van de tijdsafgeleide van het locale intensiteitsverloop 
een masker afgeleid. Waarden voor het masker worden bepaald voor afeonderlijke 
ruimtelijke posities en uit de waarde van de tijdsafgeleide van het locale intensiteitsverloop 
ter plaatse van de betreffende ruimtelijke positie in de ipdividuele dataset- Zo'n masker wordt 
toegepast om het gebied van belang uit de betreffende dataset te segmenteren. Bijvoorbeeld is 
30 het masker uitgevoerd als binaire dataset De segmentate wordt dan uitgevoerd door in de 
betreffende dataset de datawaarden over te houden bij ruimtelijke posities waarvoor het 
masker eenzelfde waarde, bijvoorbeeld de waarde 1 heeft. Het masker kan eenvoudig worden 
verkregen door een drempelfilter op de tijdsafgeleide van het locale intensiteitsverloop toe te 
passen waannee alleen de waarden van de tijdsafgeleide van het locale intensiteitsverloop 
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voor de betreffende dataset groter dan wel kleiner in vergelijking met een vooraf ingestelde 
grenswaarde worden behouden. 

In een verdere voordelige uitvoering van de uitvinding wordt voor een of meer 
van de dalasets de ruimtelijke verdeling van de tijdsafgeleide van het locale 
5 intensiteitsverloop als beeld weergegeven. Met name wordt de waarde van het m^mum van 
de tijdsafgeleide van het locale intensiteitsverloop op de betreffende ruimtelijke positie, of 
voor het betreffende blok a weergegeven. Zo'n beeld wordt wel aangeduid als 'idpslope map \ 
Dit lukt bijzonder goed door gebruik te maken van zogeheten valse kleuren waarbij de 
kleurwaarden de plaatselijke waarde van de tijdsafgeleide van het locale intensiteitsverloop 
10 vertegenwoordigen, Bijvoorkeur wordt in cardiologische toepassingen zo'n upslope map 
gemaakt voor het gebied van belang dat betrekking heeft op het myocardium. Deze upslope 
map geeft een goede weergave van de doorbloeding van bet myocardium en is daardoor een 
geschikt technisch hulpmiddel voor de cardioloog om infarctueuze gebieden en de maten van 
het infarct op te sporen. 

15 De uitvinding heeft verier betrekking op een analysesysteem om de werkwijze 

volgens de uitvinding mee uit te voerea Zo'n analysesysteem volgens de uivinding is 
gedefinieerd in Conclusie 16. Bij voorkeur is zo'n analysesysteem uitgevoerd als een 
medisch diagnostiscb werkstarion. Verder heeft de uitvinding betrekking op een 
computerprogramma voor het analyseren van opeenvolgende datasets, Zo'n 

20 computerprogramma volgens de uitvinding is gedefinieerd in Conclusie 17. Bij voorbeeld 
kan een computerprogramma volgens de uitvinding worden geladen in het geheugen van een 
analysesysteem, zoaLs een xnedisch diagnostiscb weikstation. Het computerprogramma 
volgens de uitvinding kan worden opgeslagen op een datadrager, zoals een CD-rom schijf en 
vanaf de datadrager in het geheugen van het analysesysteem geladen worden. Het is ook 

25 mogelijk om het computerprogramma vanaf een netwerk, zoals het "world-wide web ' te 
laden. 



Deze en andere aspecten van de uitvinding worden bij wijze van voorbeeld 
30 nader toegelicht aan de hand van de volgende uitvoeringsvoorbeelden en aan de hand van de 
bijgevoegde tekening waarin 

Figuur 1 een voorbeeld toont van opeenvolgende datasets waarop de 
uitvinding wordt toegepast 
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Figuur 2 toont voorbeelden van het locale intensiteitsverloop 
Figuur 3 toont voorbeelden van de tjjdsafgeleide van het locale 

intensiteitsverloop en 

Figuur 4 een voorbeeld van een 'iipslope ' map toont. 

5 

Figuur 1 toont een voorbeeld van opeenvolgende datasets waarop de 
uitvinding wordt toegepast Elk van de datasets la,lb,lc en Id in Figuur 1 stellen een 
afbeelding van een dwarsdoorsnede van het hart 10 van de te onderzoeken patient voor. Deze 
afbeeldingen zijn met name magnetische-resonantieafbeeldingen, of CT-afbeeldingen. De 
10 datasets, of afb eeldingen, in het voorbeeld van Figuur 1 zijn opeenvolgend in de tijd, dus 

afbeelding lb volgt op la, afbeelding lc volgt op lb en afbeelding Id volgt op afbeelding lc. 
Voor de eenvoud van de figuur zijn maar vier datasets getoond, maar in de praktijk kunnen 
vele tientallen of honderden opeenvolgende datasets worden gehruikt. 
In de datasets worden telkens blokken 3a3b,3c,3d, van datawaarden genomen. In Figuur 1 is 
15 in elk van de datasets steeds een van de blokken getoond, met name ligt het getoonde blok 
3a,3b,3c,3d ter plaatse van het linkerventrikel 1 1 van het hart 10. Door Gaussische filtering 
wordt voor elk van de blokken 3a^b,3c,3d de gewogen gemiddelde datawaarde berekend. 
Door deze gewogen gemiddelde datawaarden in de tijd uit te zetten ontstaat voor de blokken 
3a, 3b,3c en 3d het locale intensiteitsverloop; dit is getoond als de getrokken grafiek 12 in 
20 Figuur 2. Langs de tijd-as is hier het volgnummer uitgezet van de opeenvolgende datasets, of 
afbeeldingen, ook wel als 'frame* aangeduid. Zoals uit Figuur 2 en ook uit Figuur 3 blijkt zijn 
voor het voorbeeld opeenvolgende datasets met meer dan 60 individuele datasets of 
afbeeldingen gebruikt. Het locale intensiteitsverloop 12 dat in Figuur 2 is getoond vertoont 
een sterk maximum dat overeenkomt met de boluspassage door het linker ventrikel. De 
25 gestreepte grafiek in Figuur 2 toont het locale intensiteitsverloop 1 3ter plaatse van het reenter 
ventrikel. Hier treedt het m a ^'™"™ in de plaatselijk gemiddelde datawaarden pas bij 
ongeveer het 50ste frame op, doordat de contrastbolus later door het rechter ventrikel dan 
door het linker ventrikel passeert. 

Figuur 3 toont de tydsafgeleiden 21,22 van het locale intensiteitsverloop. De 
30 tijdsafgeleiden 21,22 hebben elk de eigenschap dat een hoge positieve maximum wordt 
gevolgd door een sterk negatief minimum. Aan de hand van deze eigenschap is het 
gemakkelijk de locale intensiteitsverlopen ter plaatse van het linker en rechter ventrikel op te 
sporen, Ter vergelijkmg is met een streep-snppellijn het locale intensiteitsverloop buiten het 
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hart getoond, hier is de helderheidswaarde, dus de datawaarde, in de afbeelding ongeveer 
constant. 

Volgens de uitviading worden voor vele blokken in elk van de data sets bet locale 
intensiteitsverloop bepaald. Door vergelijking met een referentie, Deze vergelijking wordt 
5 bijvoorbeeld uitgevoerd met gebruik van de gelijkheidsmaat Hiertoe wordt een overlap van 
de tijdsafgeleide van het betreffende locale intensiteitsverloop p met de tijdsafgeleiden van 
het reference intensiteitsverloop q bepaald: 

And > gend ) " gggfo^ i glBS > 
0 CPcnd ^ Qnmn*<ia* ^ Pvm) 

Hier geven de subsripts *end s en 'start' het gekozen eindmoment en beginmoment van het 
1 0 locale intensiteitsverloop en het referent! e intensiteitsverloop aan. De overlap O is mil 

wanneer het locale intensiteitsverloop eerder begin dan wel later ophoudt dan het referentie 
intensiteitsverloop. De overlap geeft aan in welke mate de drager Q support*) van het locale 
intensiteitsverloop in de opeenvolgende datasets, dus bijvoorbeeld in de tijd, saxnenvalt met 
de drager van het reference intensiteitsverloop met Vervolgens is een geschikte 
15 gelijkheidznaat: 

0 x£Q 

Hierin is H de Heavyside stapfunctie: H(x) = ^ x >§ >Pn**e&Pmm de waarden van het 

locale maximum en het locale minimum in de tijdsafgeleide van het locale 
intensiteitsverloop. In de toepassing in de cardiologie betreffen dit maximum en het 

20 minimum de aankomst van de contrastbolus en het verlaten van de contrastbolus. De 
parameter O m m is een minimum overlap met de tijdsafgeleide van het reference 
intensiteitsverloop. In de gelijkheidsmaat S wordt ex mee rel cen i n g gehouden dat in de 
tijdsafgeleide van het locale intensiteitsverloop het minimum na het maximum volgt. Aan de 
hand van de gelijkheidsmaat worden datawaarden in de datasets c.q. pixels in de 

25 opeenvolgende afbeeldingen geclassificeerd als behorende tot het gebied van belang, 
bijvoorbeeld de afbeelding van het linker ventrikeL 

Figuur 4 toont een voorbeeld van een 'upslope 3 map. Deze 'upslope map 1 
betreft het myocardium. Hier is vooraf het gedeeite van het myocardium dat in de 
afbeeldingen is weergegeven opgespoord als het gebied van belang. In de afbeelding van het 

3 0 myocardium is vervolgens plaatselijk de waarde van het maximum van de tijdsafgeleide van 
het locale intensiteitsverloop als 'valse* kleur of grijswaarde weergegeven- Duidelijk is 
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achtbaar dat er gedeelten 31 zijn die behoorlijk goed doorbloed zijn; daar worden behoorlijk 
hoge maximum waarden van de tijdsafgeleide van het locale intensiteitsverloop. In andere 
gebieden 32 is de doorbloeding kennelijk veel slechter, want daar worden lang niet 20 hoge 
maxima van de tijdsafgeleide bereikt. 
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CONCLUSEBS: 



1 - Werkwijze voor het analyseren van opeenvolgende datasets 

- waarin de afeonderlijke datasets dataelementen omvattexx die datawaarden toevoegen aan 
niimtelijke posities waaimee 

- een locaal intensiteitsverloop [ifr#]wordt opgemaakt uit datawaarden in opeenvolgende 
datasets in overeenkomstige ruimtelijke posities en 

- uit een of meer van de opeenvolgende datasets een gebied van belang wordt opgespoord 
aan de hand van het locale intensiteitsverloop 

- waarbij in het gebied van belang het locale intensiteitsverloop voldoet aan een vooraf 
gegeven eigenschap. 

2. Een werkwijze voor het analyseren van opeenvolgende datasets volgens 

Condusie 1, waarin het locale intensiteitsverloop wordt opgemaakt voor respectievelijke 
blokken van meerdere ( 



15 3 . Een werkwijze voor het analyseren van opeenvolgende datasets volgens 

Conclusie 1 of 2, 

- waarin het gebied van belang wordt opgespoord aan de hand van variaties in het locale 

intensiteitsverloop, in het bijzonder aan de hand van een tijdafgeleide ^( x >0 van 

dt 



locale intensiteitsverloop. 



4. Een werkwijze voor het analyseren van opeenvolgende datasets volgens 

Conclusie 3, 

- waarmee het gebied van belang wordt opgespoord door blokken van dataelementen op te 
sporen waarin 

25 - de variaties in het locale intensiteitsverloop groter zijn dan een vooraf ingestelde 

plafondwaarde en/of door blokken van dataelementen op te sporen waarin de variaties 
in het locale intensiteitsverloop kleiner zijn van een vooraf ingestelde bodemwaarde. 
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5- Een werkwjze voor het analyseren van opeenvolgende datasets volgens 

Conclusie 1, waarin 

- een classificatie van dataelementen in 66n of meer a&onderlijke datasets wordt 
uitgevoerd 

5 - waarin de classificatie van de dataelementen in de betreffende dataset(s) aangeeft of het 
betreffende dataelement al dan niet tot het gebied van belang boon. 

6. Een werkwijze voor het analyseren van opeenvolgende datasets volgens 
Conclusie 5, waarbij 

10 - bij de classificatie geschiedt aan de hand van een gelijkheidsmaal van het locale 
intensiteitsverloop met een referentie intensiteitsverloop in het gebied van belang. 

7. Een werkwijze voor het analyseren van opeenvolgende datasets volgens 
Conclude 5, waarbij 

15 - bij de classificatie geschiedt aan de hand van een correlate van het locale 

intensiteitsverloop met een gemiddeld intensiteitsverloop in het gebied van belang. 

8 . Een werkwijze voor het analyseren van opeenvolgende datasets volgens 
Conclusie 1 en waarbij 

20 - in afzonderlijke datasets gedeelten met ruimtelijke gradiSnten van datawaarden worden 
gemaskeerd voor zover in de betreffende gedeelten de niimtelijke gradiSnten ham 
modulus boven een vooraf ingestelde aanyaardbare gradientmodulus ligt 

9. Een werkwijze voor het analyseren van opeenvolgende datasets volgens 

25 Conclusie 1 waarbij de opeenvolgende datasets ten opzichte van elkaar in overeenstemming 
worden gebracht. 

10. Een werkwijze voor het analyseren van opeenvolgende datasets waarbij 

- meerdere gebieden van belang worden opgespoord aan de hand het locale 
30 intensiteitsverloop en waarin in genoemde gebieden van belang het locale 

intensiteitsverloop voldoet aan een vooraf gegeven eigenschap en 

- maviTniim intensiteitsprojecties (MIPs) worden bepaald van de respectievelijke gebieden 
van belang en 
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- uit verschUlen tussea deze maximum intensiteitsprojecties een thema-afbeelding (feature 
image) woidt gevonnd. 

1 1 . Een werkwijze voor het analyseren van opeenvolgende datasets volgens 
5 Conclusie 10, waarin 

in de thema-afbeelding een centrum van net gebied van belang wordt bepaald. 

12. Een werkwijze voor het analyseren van opeenvolgende datasets volgens 
Conclusie 1 1 , waarin 

10 - de themaafbeelding wordt getransformeerd naai poolcoordinaten het genoemd Centrum 
als oorsprong en 

- in deze getransformeerde thema een begrenzing van het gebied van belang wordt 
opgespoord. 

15 13. Een werkwijze volgens Conclusie 3, waarin 

- uit de tijdsafgeleide van het locale intensiteitsverloop voor een individuele dataset een 
masker wordt afgeleid en 

- met gebruik van het masker het gebied van belang uit de betreffende dataset wordt 



20 

14. Een werkwijze volgens Conclusie 13 waarin het masker wordt afgeleid door 

toepassen van een drempelfilter op de tijdsafgeleide van het locale intensiteitsverloop voor de 
betreffende dataset 

25 15. Een werkwijze volgens Conclusie 3, waarin vooi een individuele dataset een 

ruimtelijke verdeling de tijdsafgeleide van het locale intensiteitvexloop wordt weergegeven. 

16. Analysesysteem voor het analyseren van opeenvolgende datasets 

- waarin de afzonderlijke datasets dataelementen omvatten die datawaarden toevoegen aan 
30 ruimtelijke posities en 

- welk analysesysteem is ingericht voor het 

- opmaken van een locaal intensiteitsverloop [I(x>t)] uit datawaarden in opeenvolgende 
datasets in overeenkomstige ruimtelijke posities en 
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- opsporen uit een of meer van de opeenvolgende datasets een gebied van belang aan de 
hand van het locale intensiteitsvexloop 

- waarbij in het gebied van belang het locale intensiteitsverloop voldoet aan een vooraf 
gegeven eigenschap. 

5 

17. Computeiprogramma voor het analyseren van opeenvolgende datasets 

- waaiin de afzonderlijke datasets dataelemcnten omvatten die datawaarden toevoegen aan 
ruimtelijke posities en 

- welk computersysteem instructies omvat voor het 

10 - opmaken van een locaal intensiteitsverloop [I(x, t)] uit datawaarden in opeenvolgende 
datasets in overeenkomstige ruimtelijke posities en 

- opsporen uit een of meer van de opeenvolgende datasets een gebied van belang aan de 
hand van het locale intensiteitsverloop 

- waarbij in het gebied van belang het locale intensiteitsverloop voldoet aan een vooraf 
15 gegeven eigenschap. 
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ABSTRACT; 



The invention involves analysis of successive datasets. From these successive 



datasets a local intensity progress is formed from datavalues in successive datasets at 
corresponding positions in each of the datasets. On the basis of the local intensity progress a 
region of interest is detected in the individual datasets. In particular the time-derivative of 
5 the local intensity progress is used to detect the region of interest. The invention is notably 
useful in cardiological applications to separate the image of the myocardium from a sequence 
of 3D magnetic resonance imaging reconstruction. 



Fig.l 
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